Hawking Radiation from a Reisner-Nordstrom Domain Wall 



ON 

o 
o 

O 

Q 



Eric Greenwood 

HEPCOS, Department of Physics, SUNY at Buffalo, Buffalo, NY 1^260-1500 

We investigate the effect on the Hawking radiation given off during the time of collapse of a 
Reisner-Nordstrom domain wall. Using the functional Schrodinger formalism we are able to probe 
the time-dependent regime, which is out of the reach of the standard approximations like the Bo- 
golyubov method. We calculate the occupation number of particles for a scalar field and complex 
scalar field. We demonstrate that the particles from the scalar field are unaffected by the charge of 
the Reisner-Nordstrom domain wall, as is expected since the scalar field doesn't carry any charge, 
which would couple to the charge of the Reisner-Nordstrom domain wall. Here the situation effec- 
tively reduces to the uncharged case, a spherically symmetric domain wall. To take the charge into 
account, we consider the complex scalar field which represents charged particles and anti-particles. 
Here investigate two different cases, first the non-extremal case and second the extremal case. In 
the non-extremal case we demonstrate that when the particle (anti-particle) carries charge opposite 
to that of the domain wall, the occupation number becomes suppressed during late times of the 
collapse. Therefore the dominate occupation number is when the particle (anti-particle) carries the 
same charge as the domain wall, as expected due to the Coulomb potential carried by the domain 
walls. In the extremal case we demonstrate that as time increases the temperature of the radiation 
decreases until when the domain wall reaches the horizon and the temperature then goes to zero. 
This is in agreement with the Hawking temperature for charged black holes. 
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I. INTRODUCTION 

An important question in theoretical physics is that of 
quantum radiation from collapsing objects. An asymp- 
totic observer, watching gravitational collapse of a mas- 
sive object, will start registering radiation of quanta com- 
ing from the fields excited by the non-trivial metric of 
the background space-time. As time progresses, the ra- 
diation exhibits more and more thermal features, till fi- 
nally when the horizon is formed the radiation becomes 
completely thermal. This is in agreement with the fact 
that the radiation from a pre-existing black hole must be 
thermal, as originally shown by Hawking, see Ref.[l|. 

The typical method for describing the induced radi- 
ation of the collapsing object is to use the Bogolyubov 
method. Here, one considers that the system starts in an 
asymptotically flat metric, then the system is allowed to 
evolve to a new asymptotically flat metric. By match- 
ing the coefficients between the two asymptotically flat 
spaces at the beginning and end of the gravitational col- 
lapse, the mismatch of these two vacua gives the num- 
ber of particles produced during the collapse. Hence, 
what happens in between is beyond the scope of the Bo- 
golyubov method. Therefore the time-evolution of the 
thermodynamic properties of the collapse cannot be in- 
vestigated in the context of the Bogolyubov method. 

Recently much work has been done with exploring the 
time-dependent aspects of the induced radiation due to 
gravitational collapse, see for example Refs. [H, 0] . So far, 
however, these investigations have only included gravita- 
tional collapse of massive objects. However, as suggested 
by the second law of black hole thermodynamics, col- 
lapsing objects can contain other observable quantities, 
such as charge and angular momentum. An important 
question is, what happens when we introduce these (or a 
subset of these) additional variables into the mix? Here 



we investigate the Hawking radiation a charged massive 
object in the form of a Reisner-Nordstrom domain wall. 
Here the question is, how does the charge effect the time- 
development of the induced radiation? 

To investigate this question, we shall use the so- 
called Functional Schrodinger formalism (see for exam- 
ple Refs.0, H,01 )■ fn general, the Functional Schrodinger 
formalism yields a functional differential equation for the 
wavcfunctional, ^[g^, $, O], where g^ v is the met- 
ric, A M the position of the object, 4> is a scalar field and 
O denotes all the observer's degrees of freedom. Since 
the Functional Schrodinger formalism depends on the ob- 
server's degrees of freedom, one can introduce the "ob- 
server" time into the quantum mechanical processes in 
the form of the Schrodinger equation. The wavefunc- 
tional is then dependent on the chosen observer's time, 
hence one can view the quantum mechanical processes 
of a given system under any choice of space-time folia- 
tion. One benefit of this formalism is that one can solve 
the time-dependent wavcfunctional equation exactly us- 
ing the invariant operator method, in particular in the 
present paper. Here we study the propagation of induced 
radiation, represented as a scalar field and complex scalar 
field in the background of gravitational collapse, which is 
governed by the harmonic oscillator equation with a time- 
dependent frequency. We solve the equations of motion 
to find the time-dependent wavcfunctional for the sys- 
tem. One can then define the occupation number of the 
induced radiation using the wavcfunctional, which is the 
gaussian overlap between the initial vacuum state and 
the wavefunction at any given later time. By taking the 
occupation number to be of the form of the Planck distri- 
bution, we can then define (3 which is also time-dependent 
and fit the temperature of the collapsing object. 

In this paper we will investigate the collapse of a 
charged spherically symmetric infinitely thin domain wall 
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(representing a shell of matter). The details about the 
collapse will depend on the particular foliation of space- 
time used to study the system. In this paper we are con- 
cerned with the Cosmological view point of the induced 
radiation during the gravitational collapse, therefore, we 
adopt the view point of an outside asymptotic observer. 
Thus, we study the collapse of the domain wall from the 
view point of a stationary asymptotic observer. 



II. SET UP 



where a dot refers to the derivative with respect to the 
asymptotic observer time t. 

By integrating the equations of motion for the spherical 
domain wall, it was shown in Ref.Q that the mass is a 
constant of motion and is given by 



M = AnaR 2 



s/l + Rj - 2nGaR^ + (8) 



Q 2 



where R T = dR/dr, while a is the surface tension (energy 
density per unit area) of the wall. 



We consider a charged spherical domain wall represent- 
ing a spherical shell of collapsing matter and charge. The 
wall is described by only the radial degree of freedom, 
R(t). The metric is taken to be the solution of Einstein 
equations for a spherical domain wall with charge. The 
metric is Schwarzschild outside the wall, as follows from 
spherical symmetry 



ds 2 = - 1 - 



1 - 



2GM 
r 

2GM 
r 



r > R(t) 
(1) 



III. RADIATION, SEMI-CLASSICAL 
TREATMENT 



Here we consider the radiation given off during gravi- 
tational collapse. We will do this for two different cases. 
First, we will investigate the radiation by coupling a 
scalar field to the background of the collapsing shell. Sec- 
ond, we will investigate the radiation by coupling a com- 
plex (i.e. charged) scalar field to the background of the 
collapsing shell. 



Scalar Field 



where M is the mass of the wall, Q is the charge of the 
wall, and 



dfl 2 = d6 



2 ' sin 2 



(2) 



In the interior of the spherical domain wall, the line ele- 
ment is flat, as expected by Birkhoff's theorem, 



To investigate the radiation, we consider a scalar field 
$ in the background of the collapsing shell. The scalar 
field is decomposed into a complete set of basis functions 
denoted by {uk{r)} 



(9) 



ds 2 = -dT 2 + dr 2 + r 2 d9 2 + r 2 sin 2 



, r < R(t) (3) 



The equation of the wall is r = R(t). The interior time 
coordinate, T, is related to the asymptotic observer time 
coordinate, t, via the proper time of an observer moving 
with the shell, r. The relations are 



dT / fdR 

d7 = V 1+ IdF 



and 



dt 1 /, fdR 



where 



2GM Q 2 



R R 2 ' 
Upon taking the ratio of Eq. (j4j and Eq. 



dT _ (1 - /) ff2 



(4) 



(5) 



(0) 



one finds 



(7) 



The exact form of the functions Uk(r) will not be im- 
portant for us. We will, however, be interested in the 
wavefunction for the mode coefficients {afc}- 

The Hamiltonian for the scalar field modes is found by 
inserting Eq.©, Eq.flTJ and Eq.Q into the action 



(10) 



Since the metric inside and outside the shell have dif- 
ferent forms, we split the action into two parts 



S 



S n 



where 

Sin = 2n I dt 



R(t) 



drr 



S out = 27T / dt I drr" 

lR(t) 



f 



T{d r <S>f 



(11) 



(12) 



l-2GM/r + Q 2 /r 2 
2GM Q 2 



(dr*) 



(13) 
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where T = dT/dt is given in Eq.([7]). Using the results 
in Ref.0], one finds that in the non-extremal case (i.e. 
when the square root in Eq.® is zero) the velocity of 
the shell, in respect to the asymptotic observer, is given 
by 



fl«A/i- 



h 2 



(14) 



where h = H/4wfi, where H is the Hamiltonian of the 
system (which is a constant of motion), and 



H=l- 2-kcjGRh- 



(15) 



Using Eq.((7]) and Eq. (fT4|) . gives that T can be rewritten 
as 



T=A/1 + (1-/) 



h 2 ' 



(16) 



(17) 



If we define the horizons as 

R± = GM ± y/(GM) 2 -Q 2 



we can then see that as R — > R±, f — ► 0. From Eq.fIB]) 
we see that as R — > R±, T ~ / — > 0. Therefore the 
kinetic term in £,•„ diverges as (R — R±)~ x in this limit 
and dominates over the softer logarithmically divergent 
contribution to the kinetic term from S ou t- Similarly 
the gradient term in Si n vanishes in this limit and is 
sub-dominant compared to the contribution coming from 
S ou t- Hence the action becomes 



S ~ 2tt / dt 



1 



drr 1 1 



drr 2 (d t <5>) 2 
2GM Q 5 



(dr®) 2 (18) 



where we have changed the limits of integrations to R+ 
since we are asymptotic observers working in the regime 
R(t) ~ R+. 

Now, we use the expansion in modes in Eq. ([9]) to write 



S 



dt 



1 

"2? 



a k M kk idk< + -a fe N fcfc /afe 



(19) 



where M and N are matrices that are independent of 
R(t) and are given by 



M kk , 
N kk , 



R+ 

=4:ir I drr 2 u k (r)uk'(r) 
lo 

2GM 



(20) 



drr z 1 



u' k (r)u' k ,(r). 



(21) 



Using the standard quantization procedure, the wave- 
function tp{a k ,t) satisfies 



/in^M-^fcfc/IIfc, + ^a k N kk ,a k , 



* = i% (22) 



where 



n. 



dai 



(23) 



is the momentum operator conjugate to a k . 

So the problem of radiation from the collapsing domain 
wall is equivalent to the problem of an infinite set of cou- 
pled harmonic oscillators whose masses go to infinity with 
time. Since the matrices M and N arc symmetric and 
real (i.e. Hermitian), it is possible to do a principal axis 
transformation to simultaneously diagonalize M and N, 
see Ref. [|| . Then for a single eigenmode, the Schrodingcr 
equation takes the form 



dt 



(24) 



where m and K denote the eigenvalues of M and N, and 
b is the eigenmode. 

We re- write Eq. lpM]) in the standard form 



1 d 2 m 2/ \,2 

^W 2 + ^ )b 



tp{b, rj) 



. dipjb, rj) 
d-q 



where 



and 



dt'f 



dt' 1 



2GM 
R 
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R 2 



(25) 



(26) 



K K 



mf ml- 2GM/R + Q 2 / R 2 



1-2GM/R + Q 2 /R 2 ' 



(27) 



We have chosen to set r)(t = 0) = 0. 

To proceed further, we need to choose the background 
spacetime, i.e. the behavior of R(t). It was shown 
in Ref.0] that the classical solution for the Reisner- 
Nordstrom domain wall in the non-extremal case is given 
by 



R(t) ^R + + {Ro-R+)e 



-f-t/R+ 



(28) 



where which is given in Eq. (|17p . is the outer radius 
and 



/- = 1 - 



GM - y/{GM) 2 - 
R~+ 



(29) 



Eq.® tells us that / = 1 - 2GM/R + Q 2 / R 2 w 
f_e~f- l f R + at late times. We are mostly interested in 
the particle production during this period. 

The spacetime is static at early times and the initial 
vacuum state for the modes is the simple harmonic oscil- 
lator ground state, 



* (M = o) = (^) 



1/4 



-mujQb /2 



(30) 
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Then the exact solution to Eq. (|25|) at later times is, see 
Ref.l, 



ijj{b, rj) = e< 



>■('/) 



1/4 



cxp 



2 V P P 



(31) 

where p v denotes the derivative of ^(77) with respect to 
•q, and p is the real solution to the ordinary (though non- 
linear) auxiliary equation 



p m +uo 2 (vi)p = -3 



with initial conditions 



P(0) 



The phase a is given by 



, P„(0) = 0. 



a(r)) 



drf 



(32) 



(33) 



(34) 



Consider an observer with detectors that are designed 
to register particles of different frequencies for the free 
field <f> at early times. Such an observer will interpret the 
wavefunction of a given mode b at late times in terms of 
simple harmonic oscillator states, {<p n }, at the final fre- 
quency Co. The number of quanta in cigenmode b can be 
evaluated by decomposing the wavefunction in Eq. (j3Tjl in 
terms of the states, {tp n }, and by evaluating the occupa- 
tion number of that mode. To implement this evaluation, 
we start by writing the wavefunction for a given mode at 
time t > tf in terms of the simple harmonic oscillator 
basis at t = 0. This is given by 



where 



i>(b,t) = Y,Cn{t) V>n (b) 



db V * n (b)^(b,t) 



(35) 



(36) 



which is an overlap of a Gaussian with the simple har- 
monic oscillator basis functions. The occupation number 
at eigcnfrcqucncy Co (i.e. in the cigenmode b) by the time 
t > tf, is given by the expectation value 



(37) 



The occupation number in the eigenmode b is then given 
by (see Appendix IB"]) 



N(t,uo) 



Ujp 

Tf 



LUp' 



ujp 



(38) 



Here we will consider only the situation when (GM) 2 > 
Q 2 , since the case where Q 2 > (GM) 2 is not as impor- 
tant for the study of Hawking radiation. This is because 



in the case when Q 2 > (GM) 2 , this represents and oscil- 
lating system (due to the repulsive Coulomb potential) 
and a black hole is never formed, see Ref.0). This occurs 
since the exponential in Eq. (|28p has an imaginary term, 
which causes the oscillation of the domain wall. In the 
case (GM) 2 > Q 2 we have gravitational collapse and the 
domain wall will collapse to form a black hole. Here we 
consider when the charge is large, close to that of (GM). 
As discussed in Ref.[||, when (GM) 2 >> Q 2 , f reduces 
to B, which is given by 



B = 1 — 



2GM 



R 



(39) 



This is the case studied in Ref. @ from the view point 
of the asymptotic observer and in Ref.0] from the view 
point of the infalling observer (one who is falling in with 
the shell). 

By calculating N, it can be checked that N remains 
constant for t < and also t > tf. Hence all the particle 
production occurs for > t > tf. There is a possibility 
that the particle production is due to discontinuities in 
the derivative of R at t = 0,tf. However, as we shall 
see below, the particle number grows with increasing tf, 
while the discontinuity at t = is fixed, and that at t = 
tf gets weaker. This indicates that particle production 
occurs only during < t < tf and is a consequence of 
the gravitational collapse. Now we can take the tf — ► 00 
limit. In Appendix [S] we have shown that p remains 
finite but p v — > —00 as t > tf — ► 00, provided c^o 7^ 0. 
However, we are interested in the behavior of N for fixed 
frequency, Co. Since Co = loqc^ 1 / 2 ^ , t — > 00 also implies 
ooa — > 0. From the discussion in Appendix [5] we also 
know that p — > 00 as loq — > 0. 

Therefore the occupation number at any frequency di- 
verges in the infinite time limit when backreaction is not 
taken into account. 

In Figure [1] we plot N versus tf- / R+ for various fixed 
values of Here we chose that Q 2 = 0A 2 (GM) 2 . 

We see that the occupation number at any frequency di- 
verges in the infinite time limit when backreaction is not 
taken into account. 

We have also numerically evaluated the spectrum of 
mode occupation numbers at any finite time and show the 
results in Figure [2] for several different values of tf-/R+. 
We compare the curves in Figure [2] with the occupation 
numbers for the Planck distribution 



N P (oo) 



1 



; /3w _ 1 



(40) 



where (3 is the inverse temperature. It is clear that the 
spectrum of occupation numbers is non-thermal. In par- 
ticular there is no singularity in N at 00 — at finite 
time, there are oscillations in N. As tf-/R+ — > 00, the 
peak at uo = does diverge and the distribution becomes 
more and more thermal. 

Now from Eq.([25|). since the time derivative of the 
wavefunction on the right-hand side is with respect to 77, 
uo is the mode frequency with respect to 77 and not with 
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cJRJf- = 2, 5, 10 




t f-/R + 



FIG. 1: Here we plot N versus tf-/R+ for various fixed values 
of QR+/f- for the scalar field case, with Q 2 = 0.4 2 (GM) 2 . 
The curves are lower for higher u)R + / /_ . 



N 



tf-/R+ = 11, 12, 13 




FIG. 2: Here we plot N versus GjR+/f- for various fixed val- 
ues of tf-/R+ for the scalar field case, with Q 2 — OA 2 (GM) 2 . 
The occupation number at any frequency grows as tf-/R+ 
increases. 



respect to time t. Eq. (|26|) tells us that the frequency in 
t is 1 — 2GM/R + Q 2 /R 2 times the frequency in 77, and 
at time i/, this implies 



U) 



(/) 



f-e 



-/-*/«+ 



(41) 



where the superscript (t) on w refers to the fact that this 
frequency is with respect to t. This rescaling of the fre- 
quency implies that the temperature for the asymptotic 
observer (with time coordinate t) can be obtained by find- 
ing the "best fit temperature" and then rescaling by 
1 - 2GM/R + Q 2 /R 2 . So the temperature seen by the 
asymptotic observer is 



T = f-e-*- t/R +p- 1 (t f ). 



(42) 



(The temperature T is not to be confused with the time 
coordinate within the spherical domain wall, also denoted 
by T.) We can fit the thermal spectrum to the collapsed 
spectrum of Figure as shown in Figure [3] to obtain the 
temperature of radiation. Here we see that the inverse 



temperature decreases as tf-/R+ increase, meaning that 
the temperature increases over time. 



ln(l + l/N) 
0.08 



t /_/*+ = 15 




wR+//_ 



FIG. 3: Here we plot ln(l + 1/N) versus QR+/f- for 
tf-/R+ = 15 for the scalar field case, with Q 2 = 0.4 2 (GM) 2 . 
The dashed line shows ln(l + 1/Np) versus G)R+/f- where 
Np is a Planck distribution. The slope gives /3 _1 and the 
temperature in Eq. (|42l) . 

We thus see that in the context of the Schrodingcr 
formalism there is evidence of Hawking-likc, but non- 
thermal radiation emitted during gravitational collapse 
for the Reisner-Nordstrom domain wall before any event 
horizon is formed. 

Here we make some general comments on the results 
of this section. Comparing the results from this section 
with those found in Ref. Q, we see that the presence of 
the charge has qualitatively no effect when considering 
the radiation in the form of a neutral scalar particle. This 
is expected since the scalar field is only influenced by the 
gravity, not by the charge, since the scalar field itself has 
no charge. Therefore one would expect that there is no 
effect from the charge. There are, of course, minor quan- 
titative difference due to differences in the background 
metric (Schwarzschild versus Reisner-Nordstrom), such 
as the surface gravity. However, for this analysis they 
are not important. If one wants to investigate the ef- 
fect of the charge, one would need to consider a complex 
scalar field. This is done in the next section. 



B. Complex Scalar Field 

The above analysis is consistent with what one would 
expect when coupling a scalar field to the background of 
the Reisner-Nordstrom domain wall, since the scalar field 
has no charge. Therefore, it is advantageous to investi- 
gate the effect of adding charge to the scalar field and 
investigating the effect of the charge. This is done by 
considering a complex scalar field coupled to the back- 
ground of the Reisner-Nordstrom domain wall. Thus the 
action in Eq. (fTU|) becomes 



So, = 



(43) 



6 



where 



Dn = ft + iqAf, 



(44) 



is the covariant derivative, q is the charge of the field 
and A u is the vector potential. The complex scalar field 
is again decomposed into a complete set of basis functions 
denoted by {vk(r)} 



* = X] K(*H( r ) + ib k(t)v k (r)) 



(45) 



with all others being zero. To ensure this, we take the 
time derivative of the complex mode expansion, Eq. (|45j) 
to be 



d t <S> = -iJ2 [a k {t)v k {r) - ib k (t)v k {r) 

k 



(51) 



Now, we use the expansion in modes in Eq. (|45p and 
corresponding derivative in Eq. (j5Tj) to write Eq. (|49|) as 



Again, since the metric inside and outside the shell 
have different forms, we split the action into two parts 



S — Sir 



S 



(46) 



where Si„ and S ou t are given in Eq. (TT2|) and Eq. (fT3|) , re- 
spectively, with the replacement of the covariant deriva- 
tive in S ou t- However, from Ref.[|[, we can write the 
vector potential for the collapsing shell as 



(A J = ( ^,0,0,0 



(47) 



hence there is no radial component for the vector poten- 
tial. Thus we can write D r — > ft in Eq. (|49|) and we can 
expand Eq. lfTB"]) as 



S„ut =27T 



dt I drr 2 

R(t) 



l-2GM/r + Q 2 /r 2 



|ft*| 



[(ft*)** - (ft**)*] - ^l*| 2 



(48) 



Therefore the kinetic term in Si n diverges as (R — Rh) -1 
in this limit and dominates over the softer logarithmi- 
cally divergent contribution to the kinetic term from S ou t- 
Similarly, the potential term in 5;„ vanishes in this limit 
and is sub-dominant compared to the contribution com- 
ing from Sout- Also, if we assume that the charge of the 
induced radiation is small, we can then ignore the term 
in S ou t that is proportional to (qQ) 2 . Hence, 



S ~ 2tt / dt 



1 

7ji 



drr 2 |ft*l 



iqQ / drr [(ft*)** - (ft**)*] 
Jr + 

2GM Q 2 



drr z 1 



R+ 



|ft*| 



(49) 



where we have changed the limits of integration to i?+ 
since we are working in the regime R(t) ~ 

Here the mode coefficients satisfy the commutation re- 
lations 



[6fc(*) 1 o fc '(*)]= 6 fc (i),Mt) = s k,k 



(50) 



dt 



1 

2/ 



dkAkk'Cik' + bkAkk'h' 



qQ 



O-k^kk'^k' — bk'Bkk'bk' 



r [akCkk'a>k> + bkCkk>bk>] (52) 



where the matrices A, B, and C are matrices indepen- 
dent of R(t) and are given by 



f R + 2 
Akk' = 47r / drr VkV k > 
Jo 

Bfcfc' = 47r / drrv k v k ' 
Jo 

C k «=4n ^ drr 2 + 



(53) 
(54) 



To investigate the action of the mode coefficients, 
we consider the Noether current. The corresponding 
Noether current is given by 



7V$ =2tti 



drr 2 (ft*) - iqQ / drr* ** 
/ Jo Jr h I 



1 



drr 2 (ft*)* +iqQ / drr* * 
/ Jo Jr h J 



--2ni fn*** - n$i 



(56) 



Using the mode expansion in Eq. (|45[) and corresponding 
derivative in Eq. (|51[) we can write Eq. (|56j) as 



A$ = jA k k' (d k ak> - b k b k '^j + qQ&kk> (a k a k > + b k b k >) 

(57) 

Here we see that the the Noether current satisfies the 
commutation relations 



[iV$,o] = a, and [iV$,6] = —b. 



(58) 



Hence we see that the a modes create particles with +q 
worth of charge, while the b modes create particles with 
— q worth of charge. 

Using the standard quantization condition, the wave- 
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function ip(a k ,b k ,t) satisfies 



dip 



[^.(A-^fck-nfc- +ir k {A~ 1 ) kk ,TT kl 



and 
.dip(d,r}) 



dr) 



qQ 



[n fe B fefc »(A ^wwaw -7r fe Bfefc"(A 1 ) k » k >b k ,] 



2m <9<i 2 

i (q 2 Q 2 v 2 

2 



i — qQ 

m 



^)d 2 



k"k> 



+ C kk > a k i 



ip(d,rj) 
1 K 



- — — -t^ q Q + -- d 2 
2m dd 2 m 2 / 



1>{d,rj) (65) 



(/!^B| fc ,,(A-V* + C fcfc - 



by 







where 



it 



-i- — , and 7r fc 
da k 



db k 



(59) 



(60) 



are the momentum operators conjugate to a k and b k , 
respectively. 

To solve this we will assume that the wavefunction can 
be separated as, 



ip{c, d, t) — ip(c, t)ip(d, t). 



(61) 



where we used the fact that q is small and / — > as 
R — > R + , and where c and d are the eigenmodes of a(t) 
and b(t), and m, j/ and if are the eigenvalues of the 
matrices A, B and C, respectively. 

To solve Eqs. (|M|) and (f6"5)) . we use the ansatz for the 
c modes 

^(c, t) = el d ^Qv/ m a(c, t) = e qQm/m J:(c, t) (66) 

similarly the ansatz for the d-modes 

iP(d, t)=e~S dmQy / m a{d, t) = e -« QOT / m S(d, t) (67) 

where we used the fact that q, Q, y and m are rj indepen- 
dent. Therefore we can write Eqs. ([M|) and (|6"5|) as the 
usual harmonic oscillator equation given in by 



Separating the Schrodinger equation then gives the two 
equations 



,dcr(c, rj) 
dr] 



J_cP_ 

2m dc 2 



-w 2 (rj)c 2 



.&1>(c,t) vf 



dt 



^n^A- 1 )^^ 



i^-B kk »(A 1 )k"k'Skk> 



and 



+ ^ ( / 



q 2 Q 2 



r>2 

' D kk' 



(A + C*fc/ ay ip{c,t) 



. da(d, rj) 



2m ad 1 2 



S(c, rj) (68) 



S(d,77) (69) 



(62) 



and 

dip(d,t) 



— 7r fc (A 1 ) kk m k/ 



£^-B fefe »(A 1 ) k » k >5kk> 



where 77 and a; (77) are given by Eg. ([26]) and Eq. (f27|) . re- 
spectively. 

The spacetime is static at early times and the initial 
vacuum state for the modes is the simple harmonic oscil- 
lator ground state, 



C»- b k . 



iP(c,t) 



](c(d) , ?? = 0) = (!^) 1/4 e —o C 2 ^)/ 2 



(70) 



(63) 



where we used Eq. (|60p in the second term of each equa- 
tion. Since the two equations are only different by a neg- 
ative sign in the linear term, we can then again use the 
principle axis theorem to simultaneously diagonalize the 
matrices A, B and C. This then leads to the expressions, 
written in standard form 



Then the exact solution at later times is, see Ref.JH 

1/4 



E(c,r?) 



>(/;) 



cxp 



P 



. dip(c,r)) 
drj 



1 d 2 . y 

- 7, — t~2 + l —qQ 

Am oc z m 



1 (q 2 Q 2 y 2 

2 V 2 m 

1 d 2 



K 
1 



+ i—qQ 
2m a& m 



c 2 tp(c,rj) 
1 K 



2 / 



V>(c,7?) (64) 



(71) 

where p is given by Eq. (|32p with initial conditions given 
in Eq. (|3"3")) and the phase a is defined in Eq. (f3"4"|) , and a 
similar expression holds for the d modes. 

Note, since the total wave function is given by 
ip(c,d,t) = ip(c,t)ip(d,t), there is no effect due to the 
charge for the total wavefunction. Therefore the total oc- 
cupation number is independent of the charge, hence the 
induced radiation reduces to the uncharged case. How- 
ever, the occupation of each mode is dependent on the 
charge, due to the presence of the exponential term. 

Analogous to the uncharged case, consider an observer 
with detectors that are designed to register particles of 
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different frequencies for the free fields <j> and <f>* at early 
times. Such an observer will interpret the wavefunction 
of a given mode c and d at late times in terms of simple 
harmonic oscillator states, {tp n }, at the final frequency 
Co. The number of quanta in eigenmode c and d can be 
evaluated by decomposing the wavefunctions in Eqs. ([66|) 
and (|67p in terms of the states, {</? ra }, and by evaluating 
the occupation number of that mode. To implement this 
evaluation, we start by writing the wavefunction for a 
given mode at time t > tf in terms of the simple harmonic 
oscillator basis at t = 0. This is given by 

ip(c, d,t) = y] c n (t)c m (t)ip n (c)ip m (d) (72) 

n,m 

where 



dcifn (c]ip(c, t) and c„ 



ddtpl(d)ip(d,t) 

(73) 

which are the overlap of a Gaussian with the simple har- 
monic oscillator basis functions. The occupation number 
at eigenfrequency Co (i.e. in the eigenmode c and d) by 
the time t > tf, are given by the expectation values 



N c (t,uj) = ^n|c„| 2 and N d (t,oo) 



^ m |c m | 2 



(74) 

The occupation number in the eigenmode c is then given 
by (see Appendix [Cj 



N c (t,u>) 



2qQ m /m,UP 

V2 



1 



1 

Cop 7 



Cop 



(75) 

and the occupation number in eigenmode d is then given 

by 



N d (t,Co) = e 



= -2qQ m /m UP 



V2 



1 - 



1 

Cop 2 



Pv_ 

Cop 



. (76) 

To analyze the complex scalar field case, we wish to 
investigate two different scenarios. The first being the 
non-extremal case. In the non-extremal case there are 
two horizons, however, as far as the asymptotic observer 
is concerned the important horizon is the outside horizon. 
The second is the extremal case. In this case, the two 
horizons collapse to form one horizon and the domain 
wall can be characterized by its charge alone. 



1. Non-Extremal Case 

Here we will analyze the induced radiation in the non- 
extremal case. As discussed in the previous section, here 
we will consider only the situation when (GM) 2 > Q 2 . 

In Figures 0] and [5] we plot TV versus tf-/R+ for 
various fixed values of CoR + /f-. Here we chose that 
Q 2 = 0.4 2 (GM) 2 and q = 0.2(GM) for illustration pur- 
poses. We see that the occupation number for each mode 



at any frequency diverges in the infinite time limit when 
backreaction is not taken into account. 
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FIG. 4: Here we plot JV versus tf-/R+ for various fixed values 
of u>R+/f- for the non-extremal case, with Q 2 = OA 2 (GM) 2 
for the d modes. The curves are lower for higher uoR+ / f- . 
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FIG. 5: Here we plot iV versus tf-/R+ for various fixed values 
of uiR+/f- for the non-extremal case, with Q 2 = 0.4 2 (GM) 2 
for the c modes. The curves are lower for higher d)R+/f-. 

We have also numerically evaluated the spectrum of 
mode occupation numbers at any finite time and show 
the results in Figures [H] and [7] for several different values 
of tf- 1 R+ . We compare the curves in Figures [H] and [7| 
with the occupation numbers for the Planck distribution 
in Eq. (|40[) . where again (3 is the inverse temperature. It 
is clear that the spectrum of occupation numbers is non- 
thermal. In particular there is no singularity in N at 
oo = at finite time, there are oscillations in N. As 
tf-/R+ — ► oo, the peak at oo = does diverge and the 
distribution becomes more and more thermal. 

We can fit the thermal spectrum to the spectrum of 
Figures [5] and \7\ as shown in Figures [5] and [5] to obtain 
the temperatures of radiation 



0.21 



f_ 

R + 



0.21 



\J (GM) 2 - Q 2 



l.32T H (77) 
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tfJR+ = 11, 12, 13, 15 
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FIG. 6: Here we plot N versus uiR + /f^ for various fixed 
values of tf-/R+ for the non-extremal case, with Q 2 — 
OA 2 (GM) 2 for the d modes. The occupation number at any 
frequency grows as tf-/R+ increases. 
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FIG. 7: Here we plot N versus u)R+/f~ for various fixed 
values of tf-/R+ for the non-extremal case, with Q 2 = 
OA 2 (GM) 2 for the c modes. The occupation number at any 
frequency grows as tf-/R+ increases. 



for the c modes, and 



T « 0,15^- = 0.15 ^ {GM £ Q - = 0MT H (78) 
i£_l_ 

for the d modes, where 



ln(l + l/Ay 
0.15 

0.10 

0.051 



t /_//?+ = 13, 15, 16 




pRJf- 



1000 2000 3000 4000 



FIG. 8: Here we plot ln(l + l/A0 versus for = 16 

for the non-extremal case, with Q 2 = OA 2 {GM) 2 , for the d 
modes. The dashed line shows ln(l + 1/Np) versus u)R+/f- 
where Np is a Planck distribution. The slope gives P' 1 and 
the temperature in Eq, (|42[) . The slope of the best fit line 
decreases as tf-/R+ increases. 



0.12 



t /_//?+ = 13, 14, 16 




pRJf- 



1000 2000 3000 4000 



FIG. 9: Here we plot ln(l + 1/JV) versus QR+/f- for 
tf-/R+ = 16 for the non-extremal case, with Q 2 — 
0A 2 (GM) 2 , for the c modes. The dashed line shows ln(l + 
1/Np) versus ujR+Jf- where Np is a Planck distribution. 
The slope gives and the temperature in Eg. (|42fl . The 
slope of the best fit line decreases as tf-/R+ increases. 



H 



1 v/(GA/) 2 - Q 2 



2tt 



(79) 



where we used Eg . ([421) . We can see that the inverse 
temperature in each case decreases as increases, 
hence the temperature of the radiation increases over 
time. 

Comparing Figures [5] and [J] we see that for later times 
the occupation number is greater for the c modes as op- 
posed to the d modes. The occupation number of the 
d modes still grows as time increases, however, it is sub- 
dominant to the occupation number of the c modes. This 
is expected, since we are assuming the domain wall has a 



positive charge. Hence the positive charges (c modes) are 
repelled by the domain wall causing a greater occupation 
number for late times. 



Similarly, comparing Figures [7] and [S] we see that as 
time increased, the occupation number per frequency is 
greater for the c modes. Again, the occupation per fre- 
quency for the d modes increases as time increases, but 
it is subdominant to that of the c modes. Again this is 
expected due to the Coulomb repulsion. 
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2. Extremal Case 

Here we will analyze the induced radiation in the ex- 
tremal case. In the extremal case, Eq.© becomes 

/=(l-|) 2 . (80) 

It was shown in Ref. [|[ that the time dependence of / in 
the extremal limit is given by / f» (Q 2 /(t + Q 2 )) 2 at late 
times. 

In Figures [TOl and [Til we plot ln(l + l/N) versus ujRh 
for various fixed values of tj Rh for c modes and d modes, 
respectively. Here we compare the inverse temperature 
versus that in Figure [5] Similarly to the non-extremal 
case, here we see that the inverse temperature decreases 
as t/Ru increases. Now from Eq. ([25|) . since the time 
derivative of the wavefunction on the right-hand side is 
with respect to 77, u> is the mode frequency with respect 
to 77 and not with respect to time t. Eq. ([26|) tells us that 
the frequency in t is (Q 2 /(t + Q 2 )) 2 times the frequency 
in 77, and at time tf, this implies 

J,) ={S?)** < 8i » 

where the superscript (t) on w refers to the fact that this 
frequency is with respect to t. This rescaling of the fre- 
quency implies that the temperature for the asymptotic 
observer (with time coordinate t) can be obtained by find- 
ing the "best fit temperature" /3~ 1 and then rescaling by 
(Q 2 j(t + Q 2 )) 2 ■ So the temperature seen by the asymp- 
totic observer is 

r=(^)V(*/). (82) 

We can fit the thermal spectrum as shown in Figures 
ITUl and QT] to obtain the temperature of radiation. Here 
we see that the inverse temperature increases as t/Rn 
increase, meaning that the temperature decreases over 
time. Hence the temperature of the extremal black hole 
decreases as t/Rn increases. Once the shell crosses the 
horizon, the temperature will go to zero. This is con- 
sistent with the fact that the temperature is given by 

IV. CONCLUSION 

Here we consider the induced radiation from a col- 
lapsing Rcisner-Nordstrom domain wall using a semi- 
classical analysis. To do so we calculate the occupation 
of the induced radiation, N(t,u>). The occupation num- 
ber strongly depends on time, due to the fact that one 
measures the time dependent frequency, u>. The time 
dependence of the frequency depends on the observer's 
time, therefore one has the freedom to chose a particular 



t/R H = 800, 1100 

ln(l+l/W c ) 




FIG. 10: Here we plot ln(l + l/JV) versus QR H for t/R H = 16 
for the extremal case, with Q 2 = OA 2 (GM) 2 , for the c modes. 
The dashed line shows ln(l + 1/Np) versus ujRh where Np is 
a Planck distribution. The slope gives /T 1 . The slope of the 
best fit line decreases as t/Rn increases. 

t/R H = 800, 1100 

m+Wd) 




FIG. 11: Here we plot ln(l + l/JV) versus ujRh for t/R H = 16 
for the extremal case, with Q 2 = OA 2 (GM) 2 , for the c modes. 
The dashed line shows ln(l + 1/Np) versus loRh where Np is 
a Planck distribution. The slope gives p- 1 . The slope of the 
best fit line decreases as t/Rn increases. 

observer for the analysis. In this paper we have chosen to 
investigate the induced radiation from the view point of 
an asymptotic observer, an observer located at infinity. 

First we considered the case when a neutral scalar field 
is coupled to the background of the collapsing domain 
wall. In this case there is no qualitative difference from 
that of the uncharged (Schwarzschild) case. This is ex- 
pected since uncharged matter does not couple to the 
charge of the domain wall, hence as far as the scalar field 
is concerned the background metric is just Schwarzschild. 
There are, of course, minor quantitative differences be- 
tween the two cases, since the background metrics are dif- 
ferent. However, these differences are not important here, 
since they do not shed light on the impact of the charge 
of the domain wall on the induced radiation. Compar- 
ing with the results of Ref. 0] , we see that the results of 
the scalar field here arc consistent with the results found 
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here. 

To take the charge into account, we coupled a complex 
scalar field to the background of the Reisner-Nordstrom 
domain wall. Here we considered two different cases, first 
the non-extremal case and second the extremal case, re- 
spectively. First we will discuss the non-extremal case. 

In the non-extremal case, we find that as time increases 
the occupation number for charges with the same sign 
as the domain wall (here represented by the c modes) 
become dominant over that of the charges with sign op- 
posite to that of the domain wall (represented here by 
the d modes). This is an expected result due to the 
Coulomb repulsion, the like charges will be repelled by 
the wall causing a greater flux of these charges seen by 
the asymptotic observer, while the opposite charges will 
be attracted to the domain wall causing a smaller flux. 
By fitting the temperature, we find that the temperature 
of the induced radiation is on the same order as that of 
the Hawking radiation for late times. 

In the extremal case, we find that there is radiation 
induced due to the collapse of the domain wall, due to 
the time-dependent nature of the metric. Unlike the non- 
extremal case, the temperature of the radiation decreases 
as t/Rji increases. This is in agreement with Eq. ([75)) . 
where the black hole can obtain absolute zero. 

Here we demonstrated that as time increases, or as 
the domain wall approaches the horizon, the occupation 
number of the modes which carry charge with the same 
sign as to the domain wall become dominant over that 
of the occupation of the modes which carry charge with 
opposite sign to that of the domain wall. As discussed, 
this is expected due to the Coulomb repulsion (attrac- 
tion) due to the charges. Therefore, when the domain 
wall gets very close to the horizon, the flux of radiation 
measured by the asymptotic observer will be completely 
dominated by radiation with the same sign as that of the 
domain wall. The radiation distribution function is not 
quite thermal, though it becomes thermal when the col- 
lapsing object reaches its own horizon. We call such radi- 
ation Hawking-likc or prc-Hawking radiation (as opposed 
to thermal Hawking radiation from a pre-existing hori- 
zon). However, once the horizon is formed the Hawking 
effect takes over and the radiation is completely thermal. 
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APPENDIX A: p EQUATION 



In the range < t < tf, we do not have an analytic 
solution but we can derive certain useful properties. First 
note that in terms of rj 



(A2) 



f-(i-r}/R+y 

Then the equation for p after rescaling can be written as: 



drj' 2 



1 



(A3) 



where r)' = rj/R +1 g 
are 



f-(l-V') 9 3 
lJqp. The boundary conditions 



3(0) = 1, 



dg(0) 
<b)> 



0. 



(A4) 



The last term with the 1/g 3 becomes singular as g — > 0. 
Let us consider another equation with the 1 / g 3 replaced 
by something better behaved. For example 



d 2 h 
drf 2 



(to R+f 
with boundary conditions 



(A5) 



(A6) 



Eq. (|A5[) implies that h(rf) is monotonically decreasing 
as long as h(rf) > 0. Furthermore, it is decreasing faster 
than the solution for g as long as g < 1, since the 1/g 3 
term in Eq. (|A3j) is a larger "repulsive" force than the h 
term in Eq. l[A"5]) . So 



(A7) 



for all T]' such that g{rf) > 0. 

Eq. (|A5[) with initial conditions (|A6|) can be solved in 
terms of degenerate hypergeometric functions. For us, 
the important point is that the solution for h is positive 
for all rf and, in particular, h(l) > for all the values 
of cooR+ we have checked. Therefore g{r]') is positive, at 
least for a wide range of uqR+. 

Let gi = g(l) 7^ 0. Then the equation for g can be 
expanded near rf = 1. 



d v ' 2 



9i 



1 

71 



This shows that 



^-^B+ff ln(W) 
dr] f- 



(A8) 



(A9) 



In the range t < 0, u> is a constant and the solution to 
Eq. ([32]) is 



p(v) = 



1 



(Al) 



as rf — * 1. 

Hence p(rj = R + ) is strictly positive and finite while 
p. q {rj = R + ) = —00 for finite and non-zero ujq. Since 
g = ^/ujqp, and g — > 1 for ujq — > 0, we also see that p — > 
and p v — > as ujq — > 0. 



12 



In the range tf < t, w is a constant. However, the 
solution for p is not a constant, unlike in the range t < 0, 
since the constant solution 1 / y/u{tf) does not necessarily 
match up with p(tf—) to ensure a continuous solution. 
Yet it is easy to check that in this region N = and 
so there is no change in the occupation numbers. So we 
need only find N(tf — ,u) to determine N(t — ► oo,u). 



APPENDIX B: NUMBER OF PARTICLES 
RADIATED AS A FUNCTION OF TIME, FOR 
SCALAR FIELD 

We use the simple harmonic oscillator basis states but 
at a frequency u to keep track of the different w's in the 
calculation. To evaluate the occupation numbers at time 
t > tf, we need only set u = u(tf). So 



>(&) 



7T / 



n„{Vm^b) (Bl) 



where 7i„ are Hermite polynomials. Then Eq. (|73[) to- 
gether with Eq. (|31|) gives 



2 - 2 

ir^up z 



2 - 2 

ir^up z 



1/4 ta /• 
\/2 ,l n! 7 



1/4 e ia 



(B2) 



where 



and H n (0) = for odd n, we find the coefficients c„ for 
even values of n. 



(_l)n/2 e io 



i.j^fiL^!!. (B9) 



For odd n, c„ = 0. 

Next we find the number of particles produced. Let 



X = 



2 

1 

P 



(BIO) 



Then 



N(t,ui) = n l c " 

n— even 

2 



d \ - 

" I I v n— even 

2 d 1 



(n = 1)!! 



y/Qf?\P\ dx v 7 ! - X 2 
2 X 2 

^2 |p| (l- X 2)3/2- 

Inserting the expressions for \ and P, leads to 



N(t,u) 



up 

7f 



Pa 

up 



(Bll) 



(B12) 



APPENDIX C: NUMBER OF PARTICLES 
RADIATED AS A FUNCTION OF TIME, FOR 
COMPLEX SCALAR FIELD 



P = 1 



u \ p p 1 



(B3) 



To find /„ consider the corresponding integral over the 
generating function for the Hermite polynomials 



J(z) = / d£e- p t ' 2 e~ z +2 *« 



,-z 2 (l-2/P) 



Since 



Jdte-^H n (0 = £ JO 



Therefore 



^ = /|(i-|)" /2 ^(o). 



Since 



K n (0) = (-l)"/ 2 V2^! (n / )!! , n 



(B4) 

(B5) 
(B6) 

(B7) 
(B8) 



We use the simple harmonic oscillator basis states but 
at a frequency u to keep track of the different w's in the 
calculation. To evaluate the occupation numbers at time 
t > tf, we need only set u = u(tf). So 



</>n(c) 



man 1 / 4 e -mac 2 /2 
7T / V2"n! 



(CI) 



where 7i„ are Hermite polynomials. Then Eq. (|T3"|) to- 
gether with Eq. (|3T|) gives 
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1/4 
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where 



P=l-i ^ + ± 
u \ p p 2 



(C2) 



(C3) 



To find /„ consider the corresponding integral over the 
generating function for the Hermite polynomials 



J{z) 



d£e- p e' 2 e-* +2 * 



(C4) 
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Since 



Then 



^e" p S 2 / 2 W„(0 = -—J(z 



n=0 
d" 



dz r - 



z=0 



Therefore 



12% ( 2 

p V P 



I jl 



Wn(0). 



Since 



Wn(0) 



( _l)-/a^lClza!l > n 



(C5) 
(C6) 

(C7) 



=even (C8) 



and 7Y n (0) = for odd n, we find the coefficients c„ for 
even values of n, 
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For odd n, c n = 0. 

Next we find the number of particles produced. Let 



X = 
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P 



(CIO) 



n— even 



' rZ-v * * in 1 1 



a 2qQyri/m_ 
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r I I /v - n— even 

2 d 1 
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„2qQyr}/m_ 



^\p\ rfXv/Tf 
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Inserting the expressions for \ an d P, leads to 



N c {t,Q) = e 2 '*'/" 1 ^ 
v 2 



(C12) 

Similarly, following the same steps for the d-modes 
leads to 



^1 



(C13) 
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